Prepared for submission to JCAP 



GCG Parametrization for Growth 
Function and Current Constraints 

pu Gaveshna Gupta, Somasri Sen^ and Anjan A Sen 

Center For Theoretical Physics, Jamia Milha Islamia, New Delhi- 11 0025, India 
E-mail: gaveshna@ctp-jamia.res. in, ssen@jmi.ac.in, aasen@jmi.ac.in 



Abstract. We study the linear growth function / for large scale structures in a 
cosmological scenario where Generalised Chaplygin Gas (GCG) serves as dark en- 
ergy candidate. We parametrize the growth index parameter as a function of redshift 
and do a comparative study between the theoretical growth rate and the proposed 
parametrization. Moreover, we demonstrate that growth rates for a wide range of dark 
energy models can be modeled accurately by our proposed parametrization. Finally, 
we compile a data set consisting of 28 data points within redshift range (0.15,3.8) to 
constrain the growth rate. It includes direct growth data from various projects/surveys 
including the latest data from the Wiggle-Z measurements. It also includes data con- 
straining growth indirectly through the rms mass fluctuation as{z) inferred from Ly-a 
measurements at various redshifts. By fitting our proposed parametrization for / to 
these data, we show that growth history of large scale structures of the universe al- 
though allows a transient acceleration, one cannot distinguish it at present with an 
eternally accelerating universe. 
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1 Introduction 

The accelerating phase of the present cosmic expansion is attributed to an ideal fiuid 
dubbed Dark energy [1] [2] [3] [4] [5]. The evidences for the existence of such dark entity 
come mainly from the observations that measure the expansion rate of the universe 
H{a) = a/a at different epoch, where a is the scale factor of the universe. Some of the 
important tests which provide most precise probes for dark energy by measuring the 
expansion rate of the universe are - the luminosity distance measurement from stan- 
dard candles like Type la Supernovae [6] [7] [8] [9] [10], the measurement of angular 
diameter distance using standard rulers like horizon scale of sound at last scattering[ll] 
and baryon acoustic oscillations [12]. Yet the determination of the expansion rate alone 
cannot give the crucial information which can differentiate between different dark en- 
ergy models. That is to say, models with different physical origins but the same global 
expansion properties could not be separated. 

The other important probe is the cosmic growth which tests the evolution of 
the inhomogeneous part of the energy density. The growth of large scale structures, 
derived from the linear matter density contrast 6{z) = Sp/p in the universe, serves as 
an important companion test. This can provide significant insight into the properties 
of dark energy which could possibly remove the degeneracy between various models. 

The cosmic expansion history constrains the equation of state of dark energy. 
So, the usual practice is to parametrize the effective equation of state of the dark 
energy w{z) = ^j^^ to mimic the observational description of the expansion rate 
H[z). Similarly in the cosmic growth history, the matter density perturbation 5{z) is 
constrained. The standard approach is to parametrize the growth function / = in 
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terms of growth index 7 to mimic the evolution of the inhomogeneous energy density. 
As we will describe in the next section, this kind of proposal of parametrizing 6 in 
terms of 7 was first introduced by Peebles [13] and later very effectively expanded by 
Wang and Steinhardt [14] for dark energy models. Parametrizations of 7 have been 
done from different motivations [15]. In this article we introduce a new parametrization 
for growth index 7 motivated from Generalised Chaplygin Gas (GCG). 

GCG is a very interesting alternative to other proposals aiming to explain the 
observed accelerated expansion of the Universe [16]. It is an exotic background fluid, 
characterised by the equation of state 

_^ 

Pgcg 

where A and /3 are constants. Within the framework of Friedmann- Robert son- Walker 
cosmology, this equation of state leads to the density evolution 

Pgcg = PgcgO [As + (1 " A,)(l + zf^'^^^] ^ , (1.2) 

where Ag = with Pgcgo being the present value of pgcg- The choice of As and /3 

PgcgO 

uniquely specifies the GCG model. It is straightforward to check that As = —Wgcg{z = 
0). The parameter /3 plays an interesting role. For (1+/3) > 0, at early times the energy 
density behaves as matter while at late times it behaves as a cosmo logical constant. 
Whereas, for (1 + /3) < 0, GCG behaves as a cosmological constant to start with and 
in future it starts behaving like a matter fluid. With this, the late time acceleration 
can be a transient phenomena. All these diverse possibilities make the GCG equation 
of state a unique one. 

The GCG model has been successfully confronted with various obervational tests: 
high precision Cosmic Microwave Background Radiation data[17]. Supernova data[18], 
and gravitational lensing[19]. It was also shown that GCG can mimic dark energy 
models with phantom equation of state avoiding the pathology of violation of domi- 
nant energy condition [20]. Despite all these pleasing features, there are certain issues 
of structure formation that plagues the dark matter sector of this model[21]. Some 
alternate views are also suggested in that context [22]. However, here we concentrate 
on the dark energy aspect of the model. In recent studies, it has been shown that GCG 
parametrization for dark energy equation of state consistently performs very well with 
different model selection criteria e.g x^/dof, GoF, BIG, AIC etc, when compared to 
other parametrization for dark energy equation of state [23]. 

In this article we study the growth rate of the large scale structures in a cosmolog- 
ical scenario where GCG serves as the dark energy candidate. Following the Wang and 
Steinhardt [14] ansatz, we parametrize the growth index 7 in terms of the GCG param- 
eters [As, (3). To check the efficiency of our parametrization we compare the growth 
rate from two approaches (theoretical and parametrized) and obtain a very impressive 
agreement (less than 1% difference). Subsequently, we use this new parametrization 
to portray the growth of structures in a wide range of dark energy models and show 
that this new parametrization can successfully mimic the growth index of a variety of 
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dark energy models. We also constrain the dark energy behaviour using current data 
for growth mainly from galaxy redshift distortion and Lyman-a forest. We show that 
given our current knowledge of the present day density parameter Qmo from WMAP-7 
measurements, growth data allow both transient as well as an eternally accelerating 
universe. 

The paper is planned as follows: The study of growth rate for GCG model, the 
formulation of the parametrization of growth index and a comparison between the two 
approaches are discussed in Section 2. In Section 3, the new parametrization is used 
to depict the growth rate for various dark energy models and a comparative study is 
done to see the difference between the proposed and the theoretical one. We present 
the compiled data set of currently available observations and do a maximum likelihood 
analysis in Section 4. Finally, in Section 5, we conclude by summarising our results. 



2 Parametrization of the Growth Index 



The growth rate of large scale structures is derived from matter density perturbation 
6 = Spm./Pm in the linear regime, satisfying 



6 + 2-6- ATTGpm6 = 0. 



{2.1] 



While the equations that influence the background cosmology involving matter and 
GCG fluid are 



8nG 



{Pm ~\~ Pgcg): 



—SnGujgcgPgcg, 



(2.2) 
(2.3) 



Here, the dot is derivative with respect to time 't' and Wgcg is the equation of state for 
GCG given by. 



w„ 



(2.4) 



[^ + (l-A)(l + ^)3(i+«]- 

It is quite straightforward to change the variable in equation(2.1) from t to Ina to 
obtain 



{ln6)" + {ln6y^ + ln6' 



1 3 

2 - 2'^9cg{'^ - ^m(a)) 



:^m(a), 



where 



Pr, 



(2.5) 



(2.6) 



pm Pgcg 

With the help of energy conservation equation dp = — 3(p + p)dlna we further change 
the variable from Ina to fim(fl) and recast the above equation in terms of the logarithmic 

dlnS 
dlna 



growth factor / = jp^ as 



df 



3u^gcg^m(l ^"''''d^l 



2 ~ 2^9cgi^ - ^m(a)) 



:^mia) 



(2.7) 
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As we mention in the previous section, Wang and Steinhardt[14] had proposed an 
ansatz for the logarithmic growth factor: 



/ = (2.8) 

where 7 is termed as the growth index parameter. This is an excellent fit to / for 
various cosmological models with specific choice of 7. For flat dark energy models with 
constant equation of state Wq, •j is given by 

For ACDM reduces to ^[24]. For matter dominated open universe model 7 = |[25, 
26]. For dark energy models with slowly varying equation of state, the growth index 
parameter 7 is found to be function of Qm and w. As the growth data spans over a range 
of redshift, people have parametrized 7 as a function of redshift. In one such attempt 
[27], 7 was proposed to be parametrized as 7(2;) = 70 + 1' z (where 7' = j^{z = 0)) to 
distinguish between models of dark energy and models of modified gravity. There are 
also other proposals from different motivations [28, 29]. 

Here we parametrize 7 in terms of GCG parameters Ag and /3. For that we assume 
the original ansatz by Wang and Steinhardt[14] 

/ = fi;^(^™), (2.10) 

where the growth index parameter 7(r2m) can be Taylor expanded around Qm = 1 as 

li^m) = l\in^=i) + - 1) J^l(f^^=i) + 0{Qm. - 1)'. (2.11) 
Now, equation (2.7) can be rewritten in terms of 7 as 

(2.12) 

Once we differentiate the above equation around = 1, after some simple calcula- 
tions, we find the zeroth order term in the expansion for 7 as 

7 = 3 /^-^^-^^ , (2.13) 

(5 - GWgcg) 

which agrees with the case of dark energy models with constant Wq (equation(2.9)). 
Following the same technique, i.e, differentiating twice around Qm = 1, we find the 
first order term in the expansion of 7 



d-f . 3{1-Wgcg){l 



(2.14) 
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Substituting equations (2.13) and (2.14) in the expansion (2.11), 7, till the first order, 
can be approximated as 



3Wa 



^ i5-6wg,g) ^ ™^ 125(1-^^3 ^ ^ 



One can now put the expression for Wgcg ( given in equation (2.4)) and subsequently 
7 can be parametrized in terms of the GCG parameters Ag, P and redshift z. 

We define normalised growth function g from the numerically obtained solution 
of equation (2. 5) as 

. ||. (2.16) 

Also the corresponding approximate normalised growth function from the parametrized 
form of / (equation (2.10)) is given by 

^t.(.)=expA*^-('^)^^, (2.17) 

where 7 is parametrized by equation(2.15). 

Ansatz (2.10)along with (2.15) provide an excellent approximation to the theo- 
retical growth factor obtained by numerically solving equation(2.7). To show this, we 
plot the difference between the numerically obtained solution of equation(2.7) for the 
growth factor / of GCG model and the approximated growth factor given by equa- 
tion(2.10) and (2.15) in the left panel of Figure 1. We also show the relative percentage 
difference in normalised growth function g{z) (in the right panel) for various values of 
As and (3. The difference between the two approaches is quite impressive (< 1%) as 
can be seen in the plots. 



3 Fitting the parametrized growth index to other dark energy 
models 

In the previous section we parametrize the growth index in terms of GCG parameters 
As and /3. In this section we use this parametrized form of 7 to find the growth rates 
for different dark energy models. We match this parametrized growth rate with the 
one that is calculated directly from the model. 

3.1 Dark Energy Models with CPL parametrization 

The basic equations governing the background cosmology for Dark Energy(DE) models 



are 



a j 3 

2 



(Pm + Pde), (3.1) 



2^ + ( - ) = -MGWdePde. (3.2) 
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Figure 1. Plot of the relative percentage difference in theoretical (numerically calculated) 
and proposed growth rate as well as the normalised growth function. TOP-LEFT: In the left 

we have plotted — - for various values of Ag and /3 in GCG model. The colour codes 
are: blue dotdashed (1) for = 0.9 /3 = 0.1, green-solid (2) for A^ = 0.9 and f3 = -1.05, 
red-dashed (3) for As = 0.8 and /3 = 0.6 and black-dotted (4) for As = 0.8 and /3 = -1.02. 



TOP-RIGHT: In the right we have plotted ^ 
model. The colour codes used are same. 



— for the same values of Ag and /3 in GCG 



where we have assumed pde = WdePde- Each component of energy density satisfies the 
conservation equation 

— = — 3if(l + Wj); i = m, de. (3.3) 
Pi 

Consequently, matter energy density scales as a^^. For dark energy we use the popular 
CPL parametrization[30, 31] for equation of state (EoS) 

Wde{a) = Wo + Wa{l - a), (3.4) 

following which p^e scales as a-3(i+«"o+«'a)g-3uia(i-a)_ ^jj-^jg conveniently includes the case 
of a constant EoS with {wq = w,Wa = 0), and the ACDM model {wq = —l,Wa = 0). 
We find the linear density contrast 6 from the equation(2.1) for such a background. 

Next, we fit the approximate growth factor /, calculated from the ansatz (2.10) 
and (2.15), to the growth factor calculated for a particular dark energy model and 
find the fitting values for the GCG parameters Ag and /3. Thus a set of {As, /3) can 
represent the growth factor of a particular dark energy model. 

3.1.1 ACDM 

In the CPL parametrization, ACDM is represented by Wq = —l,Wa = 0. The growth 
factor / for ACDM is numerically calculated by solving equation(2.1) with this choice 
of wq and Wa- Although the EoS of GCG reduces to cosmological constant for the 
choice As = 1, we fit the ansatz for / given by equations (2.8) and (2.15) to the 
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Figure 2. Theoretical fittings for GCG parametrization. In the left figures we have plot- 

ted the relative percentage difference ——j: — - in order to compare the fit of the proposed 
parametrization to that of the growth factor /, that is numerically integrated from the 
growth ODE for various DE models. The figures on the right side are of the growth index 
for the corresponding models with the GCG parametrization. TOP-LEFT: Plot of the rel- 
ative percentage difference against z for various QCDM models. The lines blue-dashed(l), 
red-dotted(2), green-solid(3) and brown-dotdashed(4) correspond to the quintessence models 
given in the Table 1 respectively from top to bottom. TOP-RIGHT: 7 for various QCDM 
models for the best fit value of 0^0(0.2895). MIDDLE-LEFT: Plot of the relative percent- 
age difference against z for dark energy models with thawing scalar field following differ- 
ent potentials. Here also the lines blue-dashed(l), red-dotted(2), green-solid (3) and brown- 
dotdashed(4) correspond to the thawing models given in the Table 1 respectively from top 
to bottom. MIDDLE-RIGHT: 7 for Thawing scalar field models with different potentials 
for ftmO = 0.2895. BOTTOM-LEFT: For Dark energy models with tachyon field following 
different potentials the relative percentage difference is plotted against z. We follow the 
same representation, i.e, the lines blue-dashed(l), red-dotted(2), green-solid(3) and brown- 
dotdashed(4) correspond to the tachyon models given in the Table 1 respectively from top 
to bottom. BOTTOM-RIGHT: 7 for Tachyon models following different potentials for the 
best fit value of QmO- 
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Parameters for various QCDM models. 


(Wo, Wa) 


As 




(-0.8,0) 


0.806544 


-0.991988 


(-1,0) 


1.0084 


-1.41556 


(-0.8,-0.3) 


0.793615 


-1.83504 


(-1.1,0.22) 


1.14599 


-2.59154 


Parameters for Thawing Scalar field models. 


(A,r) 


As 




(1.0,0.0) 


0.78782 


-1.57618 


(1.0,0.5) 


0.83270 


-1.55936 


(1.0,1.0) 


0.85769 


-1.55919 


(1.0,1.5) 


0.87423 


-1.56673 


Parameters For Tachyon models 


(At,rt) 


As 




(1.0,0.0) 


0.747729 


-1.59938 


(1.0,0.5) 


0.813433 


-1.56497 


(1.0,1.0) 


0.845355 


-1.55841 


(1.0,1.5) 


0.865417 


-1.56237 



Table 1. Values for parameters Ag and /3 in GCG parametrization that fit the actual 
growth for quintessence dark energy models (QCDM) with constant and variable equation of 
state, thawing scalar field models with different potentials, and tachyon models with different 
potentials. 



numerically obtained solution and find the fitting values for As and j3 to be 1.008 and 
— 1.41 respectively. 

With this choice of A^ and (3 we find that the proposed parametrization is able 
to fit the theoretical value for ACDM to within 0.02%. The difference in growth rates 
between the two approaches is shown in Figure 2 (Red dotted line in the top left figure) 
and result of the fit is enlisted in Table 1. 

3.1.2 Dark Energy Models with constant EoS and variable EoS 

Dark energy models with constant EoS are represented with CPL parametrization as 
Wq = w and = while DE models with variable EoS would have nonzero values 
of both Wq and Wa- Going through the same procedure as we have followed above, 
we solve numerically in equation (2.7) for theoretical / for different DE models and 
fit that to the parametrized /. This gives us different fitting values for (A^,/?) for 
different DE models, which we enlist in Table 1. With those fitting values, the GCG 
parametrization fits the theoretical values very accurately ( difference is < 0.2%). The 
differences between the two /'s and the two 7's for different DE models are shown in 
Figure 2 (top panel). The fitting values are enlisted in Table 1. 
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3.2 Thawing Scalar field models 

Now we test the effectiveness of our parametrization for models where slow rolling 
scalar field or thawing scalar field plays the role of dark energy. In this model, at early 
times the scalar field is nearly frozen at w = — 1 due to the large Hubble damping. Its 
energy density is nearly constant at this stage and its contribution to the total energy 
density of the universe is also nearly negligible. But as radiation/matter rapidly dilutes 
due to expansion of the universe, the background energy density becomes comparable 
to the scalar field energy density and the field breaks away from its frozen state evolving 
slowly to the region with larger values of equation of state parameter. 
The dynamical equation for the scalar field is 

■■ ■ dV , , 

+ 3i70+ — = 0, 3.5 

where V is the potential for the scalar field and H is the Hubble parameter described 
by the background cosmology as 

?>H^ = Pm + P4>- (3.6) 

These equations describing the background dynamics can be expressed as a set of 
autonomous equations [32]: 

7; = -37^(2 - 70) + A(2 - 7<^) v^S^, (3.7) 
fi; = 3(1 -7^)1]<^(1 -fi^), (3.8) 

A; = -v^Aj(r-l)V7^- (3.9) 

(prime denotes the derivative with respect to Ina) 
where VL^, 7^, A<^ and V are defined as follows: 

It is quite straight forward to express the Hubble parameter (equation(3.6)) in 
terms of Vt^ as 

Hence, once VL^ is determined from the above set of equations for various potential V , 
we can easily find the behaviour of the Hubble parameter and thereby the evolution 
of the growth factor / for the thawing scalar field. We consider various types of 
potentials e.g V = (p, V = (j)^ , V = e'^ and V = (f)~^ , characterized by F = 0, | , 1 and 
I respectively to obtain various growth rates for the thawing model. 

To fit the parametrization we follow the same technique adopted for the dark 
energy models and consequently find a set of fitting values for (A^,/?) corresponding 
to various potentials mentioned above. The results of the fits are provided in Figure 2 
(middle panel) and Table 1. Our parametrization fits thawing scalar field models with 
different potentials convincingly ( with more than 99% accuracy). 



-9- 



3.3 Thawing Tachyon models 

Our next target is Thawing tachyon models i.e, we test our proposal against these 
types of model. The tachyon field is specified by the Dirac-Born-Infeld (DBl) type of 
action [32] 

S = j -V{(t))^l - ^^'(t)^f,(t)^/^d^x. (3.13) 
The equation of motion which follows from (3.13) is 

+ m^{l - 02) + ^(1 - 02) ^ 0, (3.14) 

where H is the Hubble parameter. In a similar fashion like the thawing scalar field 
model, the evolution equations for the tachyon can also be cast in a set of autonomous 
equations [32] 

7^ = -67t(l - 7t) + 2v/3^At(l - 7^)^, (3.15) 
n[ = ?,nt{l-^t){l-^t), (3.16) 
A; = -v^A2(l-7,)i(r-^). (3.17) 

(prime denotes the derivative with respect to Ina) 

Here, the density parameter of the tachyon field Vtt, the EoS of the tachyon field 
7t = 1 + Aj and V are defined as follows 

a = -^=. It = ^\ (3.18) 

1 + 02 

^3/2 ^0 ' " ~ " rf02 / ^0 J • 



1 dv ^ .A^"^ 

A* = -T7^:7T' r = y— -/ — . (3.19) 



As done in the thawing scalar field case once fif is solved from the above set of au- 
tonomous equations it is easy to find the behaviour of the Hubble parameter for various 
potentials. Subsequently, the evolution of the growth factor for the tachyon field rolling 
down different potentials could be solved. We consider similar potentials for tachyon 
field as well, i.e, V = (\)^ V = (\P' ^ V = e'^ and V = 0~2 ^ characterized by F = 0, |, 1 
and I respectively. 

For fitting the parametrization we follow the same path and find the set of [As, /3) 
corresponding to various potentials. Like the other cases, here also we find that our 
parametrization fits the theoretical value to within 0.55%. We enlist the results in 
Figure 2 (bottom panel) and table 1. 

We fit the growth function of various dark energy models with our fitting function 
for growth index given in equation (2.15). In figure 2, we show the difference between 
our fitting function and the actual growth index for a given dark energy model. The 
figure shows that difference is extremely small for redshifts as large as 2 = 5 which 
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Figure 3. Plot of the relative percentage difference in theoretical (numerically calculated) 

qT{-) _ f 

and proposed growth rate. LEFT: In the left we have plotted ——j' — - for GCG model with 
As = 0.9 and /3 = —1.05 for different values of ftmo- The color codes are green(solid), 
red(dashed) and black(dotted) for QrnO = 0.25, 0.28 and 0.3 respectively. RIGHT: In the 



right we have plotted 



/ 



1 



for CPL model with wq = —0.8 and Wa = —0.3 for different 



values of ^mo- The color codes used are same. 



is sufficiently high for observational data from the growth history of the universe. In 
other words, our fitting function for growth index is quite accurate for redshifts that 
can be probed by present or future growth data. 

We also study the sensitivity of the fitting procedure on cosmological parameter 
e.g Qmo- In Figure 3, we show this for GCG case as well as for the CPL parametrization 
by taking different values for Qmo- The figure shows that the fitting does not depend 
on the choice of Qmo- 

4 Observational Constraints 

In this section we constrain our model with the current observational data . We per- 
form a maximum likelihood analysis to find the best fit parameters for 7(yls, /3, l^mo)- 
The growth data is given in Table 2, along with the references. The data in Table 2, ob- 
tained from references [42] , [43] , [44] , [35] involve redshift distortion parameter /3, which 
is related to the growth rate as /3 = //6 (b is bias that measures how galaxy traces 
mass density field). The data for parameters /3 and b at various redshifts are given in 
references [34], [35], where similar analyses are performed using different parametriza- 
tions. A catch in using this data to constrain other cosmological models is that at 
various steps in the process of analysing or converting the data ACDM model is used. 
So, if one uses the data to constrain other models, in particular the modified gravity 
models like DGP, one has to redo all steps assuming that model, starting from original 
observation. The data from reference [40] does not involve f3, but is obtained from 



- 11 - 



various power spectrum amplitudes of Ly-a forest data in SDSS. We thus define 

X/(^0m,7) = 

i 

where fobs and a fobs are obtained from Table 2 while fth{zi,'y) is obtained from equa- 
tions (2.10) and (2.15). 

Table 2. Currently available data for linear growth rate fobs used in our analysis, z is 
redshift; a is the la uncertainty of the growth rate data. 



z 


fobs 


a 


Ref. 


0.15 


0.51 


0.11 


[33], [41] 


0.22 


0.60 


0.10 


[35] 


0.32 


0.654 


0.18 


[37] 


0.35 


0.70 


0.18 


[42] 


0.41 


0.70 


0.07 


[35] 


0.55 


0.75 


0.18 


[43] 


0.60 


0.73 


0.07 


[35] 


0.77 


0.91 


0.36 


[36] 


0.78 


0.70 


0.08 


[35] 


1.4 


0.90 


0.24 


[44] 


3.0 


1.46 


0.29 


[40] 



fobsjZi) - fth{Zi,l) 

Jobs 



(4.1) 



Another dependable indirect observational probe for the density contrast 6{z) is 
the redshift dependence of the rms mass fluctuation as{z). It is defined as 

a\R,z)= W\kR)A\k,z)^, (4.2) 
Jo ^ 

where 

A^{kz) = ATrk^Ps{k,z), (4.4) 

with R = 8h~^Mpc and Psik, z) the mass power spectrum at redshift z. The function 
as{z) is related to 5{z) as 

M^) = j^^Mz = 0), (4.5) 



which leads to 



stfc(^i, Z2) = —— = --— = 1 , (4.6) 
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Table 3. The currently available data for the rms fluctuation crg(z) at various redshifts and 
references. Notice that the data from Ref. [45] were obtained using the normalized by a 
growth factor 6 and in our analysis we took this into account. 



z as (T^g Ref. 



2.125 0.95 0.17 [38] 
2.72 0.92 0.17 



2.2 


0.92 


0.16 


[39] 


2.4 


0.89 


0.11 




2.6 


0.98 


0.13 




2.8 


1.02 


0.09 




3.0 


0.94 


0.08 




3.2 


0.88 


0.09 




3.4 


0.87 


0.12 




3.6 


0.95 


0.16 




3.8 


0.90 


0.17 




0.35 


0.55 


0.10 


[45] 


0.6 


0.62 


0.12 




0.8 


0.71 


0.11 




1.0 


0.69 


0.14 




1.2 


0.75 


0.14 




1.65 


0.92 


0.20 





where we have used equation(2.17). From the redshift evolution of the flux power 
spectra for Ly-a forest, one gets the measurements for as{z) [38, 39, 45]. These data 
points are given in Table 3 along with the corresponding reference sources. 
Using the data of Table 3 we calculate the corresponding xi defined as 

Xs(^0m,7) = 

i 

where 0"^^^^ , is derived by error propagation from the corresponding la errors of as{zi) 
and cr8(2;j+i) while Sth{zi,Zi+i) is defined in equation (4.6) (We follow the method 
prescribed in [24]). The combined Xtoti^om,'!) as 

Xtoti^Om, 7) = Xji^om, 7) + Xli^om, l) (4.8) 

Minimising Xtot with respect to all the parameters As,(3 and Qmo, we found the 
best fit growth rate f{z). The corresponding best fit values of the parameters are 



Sobsi^i, Zi^l) — Sth{Zi, Zi^i) 



a 



(4.7) 



^obs.i 
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Figure 4. LEFT: The cosmological data for the growth rate f{z) along with the best fit / = 
^m{z)^ and the corresponding la errors (shaded region). All the models including ACDM 
are well within the la region. Only one model with CPL parametrization has barely crossed 
the la region (green line) at low redshifts. RIGHT: Contour Plot for (As,/3) marginalising 
Xtot o'^si' ^mo- The dark grey area represents the allowed region at Icr confidence level while 
the region bounded by the solid line represents that at 2a confidence level. 



As = 0.764, P = —1.436 and Qmo = 0.2895. In Figure 4 we have presented the best 
fit for / with corresponding la errors (shaded region) along with the cosmological 
growth data of Table 2. In the same plot we have also presented all the models so 
far considered in this article with their corresponding choice of {As,P) given in Table 
1. The intention is to find the models permissible within la confidence level of the 
available data. Clearly, the figure shows that all the models including ACDM are well 
within the la region. Only one model with CPL parametrization is barely crossing the 
la region (green line) at low redshifts. 

We further obtain the constraints on our three model parameters, e.g, Ag, (3 and 
Qmo- First we marginalise over the parameter ^2^0 assuming the WMAP7 bound [11] 
on ilrno, i^mo = 0.2669 ± 0.0288), and get the best fit values for As and (3 as 0.891 
and -1.191 respectively. This can be seen from the contours presented in right panel 
of Fig.4. 

We further marginalise Xtot and /3 respectively. In figure 5 we show 

the corresponding contour plots. We marginalise Xtot o^^^ and plot the contours 
in (/3, Qmo) plane. This along with the bound from 7 year WMAP data[ll] on flmo 
results an upper limit of /3, /3 < —1.3 at 2a confidence level. As we mentioned earlier 
(1 + /3) < represents cosmological evolution with transient late time acceleration, 
hence this result shows that growth history of the universe may allow a transient 
accelerating universe. This is an interesting outcome from our analysis. 

Marginalising Xtot ^^^^ contours in {As,flrno) plane along with the bound 

from 7 year WMAP data not only allows the ACDM {Ag = 1) but also gives a lower 
bound on As, As > 0.62 at the 2a confidence level. 
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Figure 5. Contour Plot for marginalised Xtot with the bound on QmO from 7 year WMAP 
data. LEFT: Contour Plot for (/3, Omo) marginalising Xtot o^^^ ^s- RIGHT:Contour Plot for 
(Asj^mo) marginalising Xtot over /3. The dark grey area represents the allowed region at la 
confidence level while the region bounded by the solid line represents that at 2a confidence 
level. 

5 Conclusion 

It has been realised over the last few years that observations probing cosmic expansion 
alone is insufficient to identify any particular theory explaining the late time accelera- 
tion. The study of the growth of large scale structures has emerged as a complementary 
probe for differentiating between the different dark energy models. In this work we 
explore this aspect using GCG as a possible parametrization for dark energy equation 
of state. We study the growth of the matter perturbation in GCG model and then 
propose a parametrization of the growth index 7 in terms of the GCG parameters 
that would mimic the original growth history of a GCG model. The new 7 is con- 
structed following the same ansatz as prescribed by Wang-Steinhardt [14]. We test 
the proposed parametrization against other candidates for dark energy. We find that 
GCG parametrization can fit the growth evolutions of different dark energy models 
extremely well. 

We constrain our model with latest data from redshift distortion of galaxy power 
spectra and the rms mass fiuctuation (as) from Ly-a surveys. Our compiled set of 
11 data points in the growth data set (Table 2) includes the 4 latest Wiggle-Z survey 
data [35]. In the other set of 17 data points, the growth rate is derived from erg data 
from the power spectrum of Ly-a surveys. Maximum likelihood analysis with all these 
28 data for the GCG parametrization shows that the growth in various DE models 
represented by GCG parametrization falls within la allowed region. Hence, with the 
current error bars for growth measurements, it is not possible to distinguish different 
dark energy models. Moreover, marginalised analysis puts an upper and lower 
bound on the GCG parameters /3 and Ag which are respectively —1.3 and 0.62 within 
2cr confidence level. From this analysis, one concludes that the growth measurements 
allow both the transient late time acceleration as well as the concordance ACDM model. 
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As we mentioned earlier, while minimising Xtot ^^e three parameters [As, P and 

^mo) as well as minimising marginalised Xtot o^^r Qmo and Ag respectively, the best fit 
values for /3 are coming out to be less than -1. We also minimize the over Ag, P and 
Qmo in the region > 0. The Xmin this case is 6.311. But when we allow 1 + /3 < 
as well, the corresponding Xmin 6.28 ( Note that this small value of Xmin due to 
the large error bars that we have in current growth data). In terms of likelihood value, 
this difference is negligible. Hence both eternal and transient accelerating models are 
allowed by the growth data with equal likelihood and at present there is no way one 
can distinguish between the two. With future experiments having tighter constraints 
on growth history, we hope that a transient acceleration may be distinguished from an 
eternal acceleration with sufficient confidence level. 

In conclusion, we propose a GCG parametrization for the growth history of the 
Universe and tested it against variety of dark energy models within the framework of 
Einstein gravity. In future, we intend to test this parametrization in modified theories 
as well as in scalar tensor theories of gravity. 
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